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VORTEX INSTABILITIES IN 3D BOUNDARY LAYERS : 

THE RELATIONSHIP BETWEEN GORTLER AND CROSSFLOW VORTICES 


Andrew P. Bassom and Philip Hall* 

Department of Mathematics 
North Park Road 
University of Exeter 
Exeter, Devon EX4 4QE 
U.K. 

ABSTRACT 

The inviscid and viscous stability problems are addressed for a boundary layer which can 

support both Gortler and Crossflow vortices. The change in structure of Gortler vortices 

is found when the parameter representing the degree of three-dimensionality of the basic 

boundary layer flow under consideration is increased. It is shown that Crossflow vortices 

emerge naturally as this parameter is increased and ultimately become the only possible 

vortex instability of the flow. It is shown conclusively that at sufficiently large values of 

the crossflow there are no unstable Gortler vortices present in a boundary layer which, in 

the zero crossflow case, is centrifugally unstable. The results suggest that in many practical 

applications Gortler vortices cannot be a cause of transition because they are destroyed by 

the 3-D nature of the basic state. In swept wing flows the Gortler mechanism is probably not 

present for typical angles of sweep of about 20 degrees. Some discussion of the receptivity 

* 

problem for vortex instabilities in weakly 3-D boundary layers is given; it is shown that 
inviscid modes have a coupling coefficient marginally smaller than those of the fastest growing 
viscous modes discussed recently by Denier, Hall and Seddougui (1990). However the fact 
that the growth rates of the inviscid modes are the largest in most situations means that 
they are probably the most likely source of transition. 

* Research was supported by the National Aeronautics and Space Administration under NASA Contract 
No. NAS1-18605 while the author was in residence at the Institute for Computer Applications in Science 
and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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§1. Introduction, 

Recently there has been interest in the effect of boundary layer growth on instability 
mechanisms to which the boundary layers are susceptible. Here, we shall concentrate on 
the Gortler vortex mechanism which has been shown to occur in both two- and three- 
dimensional boundary layer flows over concave walls. Much of the early theoretical work 
concerned with Gortler vortices addressed the problem of the linear stability of external 
two-dimensional flows over such concave walls. Early contributions were made by, among 
others, Gortler (1940), Smith (1955) and Hammerlin (1956). Later Hall (1982a, b) argued 
that much of this early work was fundamentally flawed for all the analyses invoked the 
parallel flow approximation (which essentially assumes that the basic flow in which the vor- 
tices lie is independent of the streamwise co-ordinate and so neglects the effect of boundary 
layer growth). This approximation enables the linear stability equations to be expressed 
as ordinary differential equations but Hall illustrated that this assumption is unjustifiable 
except in the limit of small vortex wavelength, and indeed this is the explanation for the 
considerable inconsistencies in the results of the previous studies. Additionally, in the 
case of small vortex wavelength the Gortler instability may be described by an asymptotic 
structure which takes account of boundary layer'growth in a rational manner and thence 
the parallel flow assumption is rendered superfluous in the only situation in which it has 
any relevance whatsoever. 

Hall (1982a) examined boundary layer flow over the cylinder y = 0, — oo < z < oo 
where the z-axis is a generator of the cylinder and y measures the distance normal to the 
surface. The x- co-ordinate measures distance along the curved surface, which is taken to 
have variable curvature (1/b) k (y) where b and Z are lengthscales. The Reynolds number 
Re, the curvature parameter 6 and the Gortler number G are defined by 

Re=—, 6 = 1 G = 2Re>6, (1.1) 

v b 

where Uq is a typical flow velocity in the x direction and v the kinematic viscosity of 
the fluid. Hall (1982a) investigated the flow characteristics when Re is large and 6 is 
small such that in the limit <5 — ► 0 G is held fixed at an 0(1) value. By scaling the 
spanwise co-ordinate z on the boundary layer thickness, it was demonstrated that for 
a non-dimensional vortex wavenumber 1 linearised vortex modes within a two- 

dimensional basic flow are neutrally stable at a Gortler number G = fifo e_4 + • • •> where go 
is a known 0(1) constant whose precise value is dependent upon the properties of the basic 
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flow under consideration. For larger wavelengths the problem is fully nonparallel and the 
linear stability equations, which now take the form of a set of partial differential equations, 
have to be solved numerically, see Hall (1983). This paper showed two significant features 
of this nonparallel flow problem, namely that the ideas of a unique stability curve and 
of unique growth rates at a specified downstream location are inapplicable to the Gortler 
problem because the location where a vortex commences to grow is dependent upon the 
position and the shape of the the imposed disturbance. 

Questions concerning the development of nonlinear nonparallel vortices within growing 
boundary layers were addressed by Hall (1988). This numerical investigation showed that 
as the nonlinear disturbance evolves the perturbation energy becomes concentrated in the 
fundamental and mean flow correction; a conclusion consistent with the weakly nonlinear 
theory of Hall (19826) valid for small wavelength vortices. It is well known that Gortler 
vortices set up in an experiment conserve their wavelength as they move downstream. 
Since the boundary layer itself thickens it follows that the local nondimensional vortex 
wavenumber becomes large as the vortex develops. Thus the small wavelength limit in 
the external Gortler problem is appropriate to the ultimate development of any fixed 
wavelength vortex and hence sufficiently far downstream in many flows the asymptotic 
work of Hall (1982a, 6) becomes applicable. 

As with all weakly nonlinear investigations, the results of Hall (19826) are valid only 
within a neighbourhood of the point where the imposed perturbation is neutrally stable. 
For vortices of wavenumber e~ 1 1, their development downstream of the point of neutral 

stability is governed by the solution of a pair of coupled nonlinear partial differential equa- 
tions which adopt a simple asymptotic structure at large values of X, where eX ( X = 0(1)) 
denotes the distance of the vortex downstream of the neutral point. Formally, for large 
X, Hall &: Lakin (1988) showed that this asymptotic structure could be used to deduce 
the flow configuration for fully nonlinear Gortler vortices, at which point the mean flow 
correction generated by the presence of the vortices is as large of the basic (undisturbed) 
flow itself. The vortex structure derived by Hall & Lakin (1988) essentially consists of a 
core region in which the vortex is concentrated and which is bounded by two thin layers 
in which the vortex activity is reduced to zero exponentially. Further work by Hall & 
Seddougui (1989) has shown that these thin layers are susceptible to secondary instabili- 
ties which take the form of travelling waves confined within these layers. This theoretical 
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investigation, together with the weakly nonlinear account of this form of secondary insta- 
bility described by Seddougui &; Bassom (1990), provides good qualitive agreement with 
several experimental observations, notably those of Peerhossaini & Wesfreid (1988a, b). 

All the work described above has addressed problems which arise when Gortler vortices 
occur in two-dimensional boundary layers but in many practical situations in which Gortler 
vortices are known to arise the basic boundary layer is three-dimensional. For example, 
in the case of a boundary layer flow over a three-dimensional obstacle or the flow over a 
turbine blade the three-dimensionality of the basic flow is potentially crucial and should 
not be neglected. Most significantly, the development of laminar flow airfoils has given 
rise to designs which have two areas of concave curvature on the lower side of the airfoil 
and when the wing is swept the boundary layer flow is fully three-dimensional and the 
previously mentioned analyses are largely inapplicable. 

The first attempt to describe this three-dimensionality effect was given by Hall (1985) 
who examined the Gortler mechanism in flow over an infinitely long swept cylinder. The 
results obtained were quite general and did not require a precise description of the particu- 
lar boundary layer under investigation. Hall (1985) showed that it is the relative size of the 
crossflow and chordwise flow over the cylinder which is critical in determining the vortex 
structure. He demonstrated that as this ratio, say A, was increased from zero the first 
significant change in the vortex structure from that in the two-dimensional case occurs 
when A ~ Re~ » , where Re is the (large) Reynolds number defined in (1.1). Then the vor- 
tices become time dependent and unlike the two-dimensional case, the high wavenumber 
modes no longer have vortex boundaries aligned with the flow direction. Indeed, as the 
crossflow increases further, the neutral vortices have axes perpendicular to the vortex lines 
of the basic flow. The neutral Gortler number for the vortices was predicted by a large 
wavenumber asymptotic analysis, the results of which suggested that for 0(1) values of 
the ratio of the crossflow and chordwise velocity fields the Gortler mechanism is probably 
unimportant compared with Tollmien-Schlichting and crossflow type instabilities. 

An investigation into the effect of crossflow on weakly nonlinear Gortler vortices was 
made by Bassom (1989). Rather than concentrate on vortices within a growing bound- 
ary layer with the inherent difficulties of providing an adequate description of boundary 
layer growth, that work considered vortices within a curved channel. This has the ad- 
vantage that the problem can be described by solution of ordinary differential equations 
as opposed to partial differential equations. Bassom (1989) showed that for this problem 
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equilibrium, streamwise independent weakly nonlinear vortices could only persist below a 
certain threshold value of the crossflow. 

One aim in the present paper is to extend some recent work by Denier, Hall & Sed- 
dougui (1990) (hereafter referred to as DHS) into a three-dimensional setting. The prin- 
cipal aim of DHS was to investigate the problem of providing a rational theory for the 
receptivity of Gortler vortices and specifically how the vortices may be triggered by wall 
roughness elements. In the course of this work, DHS reconsidered the stability of a vor- 
tex in a high Gortler number flow by implementing a linear, spatial stability analysis. In 
particular, for G 1, they considered the structure of vortices in two different wavenum- 
ber regimes: they reworked the analysis of Hall (1982a) for high wavenumbers 

and also examined 0(1) wavenumber vortices. This latter disturbance mode is governed 
by inviscid equations. By considering the region between these two wavenumber regimes, 
DHS identified a new structure which is relevant for vortices of wavenumber o(g») and 

which has the property that the vortices are trapped in a thin layer of thickness O (g - * ) 
at the wall. This is in contrast to the structure identified previously by Hall (1982a) who 
demonstrated that vortices of wavenumber o(g<) are confined to a thin layer embedded 
within the boundary layer at some well-defined position. 

Of paramount importance, DHS showed that within the o(g») wavenumber regime 
there exists a unique most unstable Gortler vortex according to linear stability theory. 
Hence the most dangerous mode within a two-dimensional boundary layer has been found 
and our objective now is to extend the work of DHS to consider the effect of introducing 
crossflow into the boundary layer flow. 

As previously mentioned, DHS conducted a purely spatial analysis of the instability 
modes. Here, we firstly add a crossflow into the situation considered by DHS and assess its 
effect; we continue to use the spatial approach. For the 0(1) wavenumber modes at G 1 
we show that crossflow first has a significant effect on the two-dimensional results once 
this parameter becomes 0(Re~ iG>). By examining the properties of the vortices when 
subject to this size of crossflow we can identify several limits of note. We show that as the 
crossflow increases the vortex structure takes on an identity which is essentially that of a 
crossflow instability, a mechanism first investigated by Gregory, Stuart &: Walker (1955). 
Additionally, if both the crossflow and the wavenumber of this inviscid vortex increase 
appropriately we can obtain a mode whose structure is dominated by viscous effects. It is 
found that when this occurs the ratio of crossflow to chordwise flow becomes 0 (Re~ * G* ) 
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and the (large) wavenumber of the vortices is that considered by DHS and which pertains 
to the most unstable Gortler vortex in a two-dimensional flow. Consequently, we have 
the basis for explaining how the stability characteristics of the modes found by DHS are 
changed by the presence of crossflow. We demonstrate that the introduction of crossflow 
into the problem can have a stabilising effect, at least according to our linear theory. In 
particular, whereas in DHS it was shown that stationary vortices are necessarily unstable 
at o(g») wavenumbers, this is no longer true once crossflow terms are introduced. In 
addition, for certain crossflow values there exist neutrally stable vortex modes whilst at 
large enough values of the crossflow no vortex instabilities induced by centrifugal effects 
persist. 

We then develop our description of the disturbance motion by allowing the vortices 
to be time dependent. Unfortunately, as in DHS, for all the problems treated here the 
equations which determine the vortices form a pair of coupled ordinary differential equa- 
tions with associated boundary conditions and this system is only amenable to a numerical 
solution for 0(1) parameter values. We study this numerical problem for a number of par- 
ticular parameter choices and also investigate the asymptotic description of the vortex 
structure in both the limits of large and of small vortex wavelengths. 

The procedure for the remainder of the paper is as follows. In the coming section 
we formulate the problem at hand and obtain the differential equations which determine 
the (essentially inviscid) structure for vortices of 0(1) wavenumber. After considering the 
solutions to this problem, we move op in section 3 to examine the viscous structure which 
describes the large wavenumber 0(G» ) vortices in the presence of crossflow. At this stage 
we include both the crossflow and temporal terms described in the preceding paragraphs 
and then subsequently we study particular cases. In section 4 we briefly consider the zero 
frequency crossflow problem and in the ensuing section examine the effect of allowing the 
disturbance to have periodic temporal dependence. We then proceed to investigate various 
asymptotic limits which describe the behaviour of the vortex modes in both low and high 
wavenumber cases and conclude with some discussion. 


§2. Formulation. 

Our aim here is to obtain the equations which determine the structure of high wavenum- 
ber vortices in a slightly three-dimensional boundary layer. As in Hall (1985), we consider 
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a boundary layer flowing over the cylinder y = 0, — oo < z < oo so that the z-axis is a 
generator of the cylinder and y measures the distance normal to the surface. The x co- 
ordinate measures distance along the curved surface, which is supposed to have variable 
curvature (1 jb)\ (7) where b and l are length scales. We define the Reynolds number Re , 
the Gortler number G and the curvature parameter 6 as in (1.1) where Uo is a typical flow 
velocity in the streamwise direction. The Reynolds number is assumed to be large, whilst 
6 is sufficiently small so that as 6 — ► 0 the parameter G is fixed and is order one. Again 
following Hall (1985), the basic three-dimensional boundary layer is taken to be of the 
form 

u = Uo[u(X t Y),Re-H(X,Y) > Re-?\'w(X,Y)} (l + O (.Re" *)) , (2.1a) 

where 

X = Y=(*)Re~±, (2.1 b) 

and A* is supposed to be 0(1). To compute this boundary layer profile in practice, U and d 
are determined by numerically integrating the two-dimensional boundary layer equations 
and the crossflow w is then found from the span wise momentum equation, 

We now define the variable Z = Re*z/l and perturb the basic flow. If t is a time 
variable scaled on I/Uq and if E = exp (iaZ) then we disturb (2.1a) by writing 


u = Uo [u + 6U(t,X,Y)E,vRe~^ + 6Re-±V(t,X,Y)E, 

\*Re~tiD + 6Re->W(t,X,Y)E\ (l +0 (fle"*)) , 


( 2 . 2 a) 


where 6 < 1. Similarly, we perturb the basic pressure p by writing 

p = p + Re~ 1 P(t, X, Y)E. (2.2 b) 


Substituting the expressions (2.2) into the governing continuity and Navier-Stokes equa- 


tions and linearising yields the perturbation equations 

U x +V y + iaW = 0, (2.3a) 

U t + uUx + U Ux + 'dUy + V Uy + A* wiaU = Uyy — , (2.36) 

v t + ttV x + UVx + flVy- + V Vy + A* ffiiaV + Gx(X)uU = — Py + Vyy — a 2 V, (2.3c) 
W i + u W x + A *UHix + + \*Vi d Y + A' wiaW = —iaP + W Y y — a^W. (2.3d) 
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The sc ali ngs (2.2) are those appropriate to Gortler vortices (see Hall (1982a) for 
futher discussion). We note the key features that the streamwise velocity disturbance is 
asymptotically larger than the other components and that with these scalings the pressure 
term drops out of the streamwise momentum equation (2.36). Further, viscous effects 
are brought into the perturbation equations (2.3) by the choice of spanwise and normal 
scalings. 

At this stage we consider the spatial inviscid Gortler modes of 0(1) wavenumber; in 
the context of the viscous equations (2.3) this corresponds to taking the limit G 1. DHS 
showed that for A* = 0, a = 0(1) the vortex structure is found by considering expansions 
of the forms 

U = {ito(x,F) + G~ » ui(x,Y) + . . .} exp (g> J fidxj , (2.4) 


together with similar expansions for G *V,G T W and G' 
(2.3) yields the problem for the amplification rate 0 


vo 


Uyy 


•u 


-VO 


— a 


with associated boundary conditions 


1 P. Substitution of (2.4) in 


XU y vq, (2.5 a) 


uo = 0 at Y = 0, vo — ♦ 0 as Y — > oo. (2.56) 

DHS observed that the most unstable spatial mode which is a solution of (2.5) may 
be written in the form 

vq = U exp ( — aY) , p 2 = ~-, (2.5c) 

which is valid for all a. Thus we have the somewhat remarkable result that, so long as U 
vanishes at the wall, the fastest growing inviscid mode has growth rate independent of the 
basic state. 

We are concerned with inclusion of sufficient crossflow in the problem so that the 
governing equations (2.5) are altered by this effect. We find that this occurs when A* = 
O(Ga) and then we write 

(U, V, W, P) = { U 0 (Y) + . . . ,G» Vo(Y) + ..., G±Wo{Y) + ..., GP 0 (F) + . . .} 

x exp (g» J fidx'j , 

(2.6a) 
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A* = AG>. 


( 2 . 66 ) 


At this stage we concentrate on a purely spatial approach so that on substituting (2.6) 
into (2.3) we obtain, in turn, 


(PU 0 + iaW 0 ) + ^ = 0 , 
i^(3u + iXuiv^ Uq + Sy Vo = 0, 
(fin + ilaio) Vo + x®Uo = 

(/?0 + iXatt) Wq + A tZ)y Vo = — iaFo- 


(2.7a) 

(2.76) 

(2.7c) 

(2.7<f) 


^From (2.7) we can easily show that 

— (— ipu + law) (— i j8%y + Xa®yy ) Vo = a 2 xQQy Vo- 

(2.8a) 

The above equation then controls the inviscid growth of vortices in a weakly three- 
dimensional boundary layer; it is of course the appropriate generalisation due to wall 
curvature of the well-known Rayleigh equation for unidirectional flows. 

Equation (2.8a) needs to be solved subject to suitable boundary conditions. Firstly, 
we demand the typical inviscid condition at the surface of the cylinder y — 0 so that 


( -ip a + law) ( - a 2 V 0 ) 


V 0 = 0 at Y = 0. 


( 2 . 86 ) 


In addition, to ensure that the disturbance is confined within the boundary layer, 

Vo — > 0 as Y — > oo. (2.8c) 


These boundary conditions, when added to (2.8a), constitute an eigenproblem for the 
spatial growth rate p in terms of the scaled wavenumber a and the scaled crossflow A. 
We notice that since this problem is a localised one about some point r = i, we may 
scale the curvature parameter x out of th e problem by appropriately redefining p and X. 
Consequently, for the remainder of the work described here we shall take x = 1. 

To fix our ideas upon some definite boundary layer flow we chose to consider (2.8) 
when u is given by the Falkner-Skan profile 

n = f'(Y ): /'" + //" + | (l - (/') 2 ) =«, m = f\ 0) = 0, /(oo) = 1. (2.9a) 
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This flow is closely related to the steady flow in the boundary layer along a surface of 
revolution near a forward stagnation point, see Jones & Watson (1963). Cooke (1950) 
showed that for a flow over an infinite yawed wedge at zero angle of attack then if the 
streamwise basic flow is given by (2.9a) then the spanwise component w satisfies 

w = g: g + fg = 0, g( 0) = 0, g(oo) = 1. (2.96) 

The dependences of U and XD on Y are depicted in Figure (2.1). We observe that the 
Blasius profile would not be a suitable choice for 0 and XD for then we retrieve the scenario 
described in Hall (1985). Hall showed that in flows with zero pressure gradient and with U 
and XD therefore linearly related the whole problem of Gortler vortices in three-dimensional 
boundary layers becomes degenerate and can be reduced to a two-dimensional case. Con- 
sequently, we needed to focus our attention on other physically realistic three-dimensional 
flows and the choice (2.9) was motivated by these requirements. 

Given this basic flow we then solved (2.8) for a variety of parameter regimes. We 
firstly considered the effect of increasing the crossflow A on the vortex growth rate j3 r (0 r 
denotes the real part of the complex number /?). .The results are illustrated in Figure (2.2) 
for a selection of A values and we can make some immediate deductions. When A = 0 we 
of course return to the case of DHS and /? is given by (2.5c). However, as A is increased 
then for most prescribed wavenumbers a the corresponding growth rate decreases so that 
the effect of increasing the crossflow is to stabilise the vortex flow. 

For small crossflows A (less than some critical value, say A e « 4.7) we observe that 
although the vortex is more stable than in the two-dimensional base flow case, we still 
have 0 T > 0 so that the vortex remains unstable for all a. Once A reaches the value A c 
however, Figure (2.2) indicates that there exists a region of wavenumber space in which 
vortices can no longer persist; indeed, for A > A c there exist cut-off values ai(A) and 02 (A) 
such that the vortices only exist for a < a\ or a > 02 - As A increases yet further, we 
see that ai decreases and 02 increases and in addition the growth rates p r in the lower 
wavenumber regime grow. In Figure (2.3) we have illustrated the dependence of 0 T upon 
A for those values of a less than ai(A). We can infer from this figure that as A — ► 00 the 
value ai(A) converges to some non-zero value, say a , and this behaviour is explained in 
the asymptotic work elucidated shortly. 

Figure (2.4) shows the behaviour of the imaginary part of /?, /?<, as the wavenumber a 
and the crossflow A vary. We immediately deduce from Figure (2.4a) that, for 0 < A < A c , 
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/?, is negative for all a and becomes increasingly more negative as a — ► oo. In addition, we 
know from the solution (2.5c) of DHS that /?< = 0 for all a when we consider a purely two- 
dimensional flow (X = 0). For values of X > A e , see Figure (2.46), we find, as previously 
mentioned, that the unstable vortex cannot persist in some wavenumber range which 
increases in extent as A grows. In Figure (2.46) this feature is indicated by the gap in 
the possible values of /?< for X > X c . As X is increased yet further, see Figure (2.4c), 
becomes progressively more negative for any prescribed wavenumber a and here we have 
plotted the dependence Pi(a) for a < ai(A). 

We turn now to consider various asymptotic solutions of (2.8). Firstly, we study the 
case A » 1, a < ai; i.e. the left-hand solution branches in Figures (2.26), (2.46) and those 
branches sketched in Figures (2.3), (2.4c). If we write 

P = (3X + o(X), (2.10a) 

and substitute this expansion into (2.8a) we obtain at leading order 

(— ifin + am) (^j 73 T — a2 ^o) — {—ipu + aw ) Vo = 0. (2.106) 

Since curvature efects are now negligible at zeroth order this is a standard Rayleigh equa- 
tion and we now have a connection between our three-dimensional vortex flow and the 
work of Gregory, Stuart & Walker (1955) (hereafter referred to as GSW) who were con- 
cerned with a description of stationary crossflow vortices formed in the boundary layer 
above a rotating disc. This instability is due to the inflectional character of the basic 
velocity profile and GSW showed with china-clay techniques that the crossflow instability 
took the form of a regularly spaced pattern of spiral vortices which was stationary relative 
to the disc. Stuart (in GSW), using inviscid theory, suggested that the instability could be 
associated with a particular inflectional profile in which the inflection point coincided with 
a point of zero velocity somewhere in the flow. His calculation gave the predicted number 
of vortices to be approximately four times the observed value of about 30, but the angle 
of 13° between the axes of the vortices and the radius vector was in excellent agreement 
with their experiments. 

We can appeal to Stuart’s inviscid analysis and note that there is a solution of (2.106) 
with boundary conditions (2.86, c) for which a is real and /? purely imaginary if there exists 
a point Y ^ 0 at which 

(-ipn(Y) + aw(Y)) = (- i(3u"(Y ) + aw" (Y)) = 0 . ( 2 . 11 ) 
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Some numerical work shows that for the basic profiles (2.9) a suitable solution of (2.11) 
may be found and then ft /a » -0.85. Consequently, we have, recalling (2.10), neutrally 
stable crossflow vortex type solutions of (2.8) if A > 1, a = 0(1) and 

4 ~ -0.85. (2.12) 

a A 

On figure (2.4c) we indicate, for each crossflow parameter considered, those points for which 
(2.12) is satisfied. We can conclude that the asymptotic result (2.12) is very accurate for 
surprisingly modest values of A. Once conditions (2.11) & (2.12) are satisified, the Rayleigh 
equation (2.106) has regular singular points at Y = 0 and Y = Y and the solution of 
this equation with boundary conditions (2.86) forms an eigenvalue problem for the GSW 
wavenumber a. We recall that earlier we observed that for A > l c the vortex mode can only 
persist for wavenumbers a less than ai(A) and greater than 02 (A) and that as A * 00 , ai 
tends to a non-zero constant ai m . It is now easy to see that ai^ is the GSW wavenumber 
a. We have chosen not to compute the precise value of a here for its value is dependent on 
the choice of the basic boundary layer flow and our present purpose is to concentrate on the 
underlying description of vortices in three-dimensional flows as opposed to concentrating 
on profiles (2.9) in particular. However, from Figure (2.3), we can deduce that for the 

basic flow (2.9), a « 1.3. 

Thus within the context of inviscid stability theory we now see the relationship between 
Gortler and crossflow vortices. For a two-dimensional boundary layer which is centrifugally 
unstable the growth rate of the most dangerous mode is given by (2.5c), the exact solution 
found by DHS. When the crossflow is increased from zero a maximum value of the growth 
rate develops at some wavenumber and the growth rate curve to the right of this point 
eventually crosses the zero growth rate axis. Thus at a finite value of the crossflow a 
finite band of unstable wavenumbers persists to the right of a = 0. At this stage we can 
identify this unstable band of modes as crossflow vortices since at even larger crossflows 
they reduce to the modes of GSW. At larger values of the wavenumber an unstable band 
of modes persists up to a = 00 for all values of the crossflow; it is perhaps appropriate 
to think of these modes as the remnants of the Gortler mechanism. However we shall see 
below that at high wavenumbers viscous effects become important and that when these 
effects are allowed for these remnants vanish completely at large enough crossflows. 

Having examined the asymptotic description of modes for which A > 1, a = 0(1) we 
next consider vortex states at high wavenumbers a. In Figure (2.5) we show the dependence 
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of the eigenvalue (3 of (2.8a) for a range of crossflows in the cases a = 5 and a = 10. We 
observe, as already commented upon, that increasing A tends to decrease the growth rate 
fl r and we expect that for each fixed a > 1 there exists a corresponding crossflow value 
at which the vortex is completely stabilised by the crossflow. This deduction is consistent 
with the results previously illustrated in Figure (2.2) where as the crossflow increases the 
wavenumber range over which vortices are unable to persist steadily widens. We can also 
conclude from earlier results that for a » 1 and X not large the vortex flow will remain 
unstable for as a — ► oo^, p T ~ a* when A = 0 and it is not expected that an 0(1) increase 
m the scaled crossflow A could result in a reduction in the growth rate by this dramatic 
amount. Consequently, we anticipate that near neutral vortices will be described for a > 1, 

A > 1. 

When a > 1, following DHS it is deduced from (2.8) that the vortices will be confined 
to a thin region of thickness 0(a *•) adjacent to the cylinder wall y = 0. We therefore find 
it convenient to write 

£ = aF - (2.13a) 

and suppose that around y = 0 the basic flow velocities have regular Taylor expansions of 
the forms 

u = ^na -1 £ + ~Ml 2 a~ 2 { 2 + . .., w = y 2 ia~ ^ + \y 22 a~ 2 £ 2 + . . . . (2.136) 

When A is gradually increased we find that the stability characteristics of the vortices 
are first altered significantly once A ~ 0(a>). If we write A = a»A, A = 0(1) we are led 
to the expansion 

p = i (a»/? 0 + a* /?i + ...) , (2.13c) 

with p 0 and Pi real. If in turn we expand F 0 = V W + a^Voi + ... and substitute into 
(2.8a) we obtain at the leading two orders the relationships 


PofJ-U + A/i21 = 0, 


(2.14a) 


and 


[/?imi£ + i (fioft 12 + A//22) £ 2 ] ( - Foo) - 

1 , I (2.146) 

PlVlli + - \PofJf 12 + \H 22 ) q [Pofiu + A/Z22J Foo = (mii) 2 ^Foo, 


12 



which needs to be solved subject to Voo - 0 at f - 0 and as £ * 00 (recall that we scaled 

X = 1). This problem needs a fully numerical solution but here we merely remark that 
this equation has a critical layer structure at the position 

£ _ 2 /?iMH (2.14 c) 

{PlVl 2 + ^ 22 ) 


and that this structure has to be fully analysed in order that the numerical scheme be 
effected. We are presently considering this problem in further work and hope to be able 
to report on this in the near future. However the structure outlined above is certainly 
consistent with our numerical calculations which suggest the emergence of a critical layer 

structure when a approaches the cut-off wavenumber a v 

For the present purposes, we can infer that Gortler vortex states require, at wavenum- 
bers a > 1, a crossflow of size 0(a±) in order to have a significant effect on their stability. 
Further, in this case the whole vortex structure is compressed into a thin 0{a- x ) region 
close to the cylinder surface. Of course, in the development of the inviscid equations (2.7) 
which describe the vortex we found that for 0(1) wavenumber vortices in a crossflow of 
size O(G’) (see (2.6ft)) the viscous terms are negligible. However, as a ► 00 the extent 
of the structure reduces so that at some stage we expect that the viscous terms become 
significant and no longer negligible. Guided by the asymptotic results (2.13) & (2.14), 
study of the fundamental equations (2.3) suggest that viscous terms are crucial and so a 

new regime is acheived when 

a=0(G>), (2.15a) 


and then ( 3 v 

A = o(g^) and (3 = O \G 10 ) . 


(2.15ft, c) 


Recalling the original scalings (2.4), (2.6) we then have a viscous structure in operation for 
0(G» ) wavenumber vortices in crossflows of size 0(G» ) and then the vortices develop on a 
streamwise lengthscale of 0(G~ *). Interestingly, DHS identified the unique most unstable 
vortex within a two-dimensional base flow as occuring within this wavenumber band and 
developing on this lengthscale. 

In this section we have restricted ourselves to calculations pertinent to stationary 
vortices. Much of the above analysis can be modified in a straightforward manner to ac- 
count for time-dependence of the vortices. However, our objective has been to concentrate 
on understanding the fundamental connections between the inviscid vortex modes and the 
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GSW crossflow instability and viscous, high-wavenumber vortex structures. There are nat- 
ural non-stationary counterparts of the GSW stationary instability (see Bassom & Gajjar 
(1988) for an example) and these will match with time— dependent inviscid vortices in a 
manner virtually identical to that described here. Further remarks relating our inviscid 
modes to the crossflow instability mechanism will be made in the conclusions. However, we 
now move on to examine the salient properties of the high wavenumber viscous structures 
suggested by scalings (2.15). 
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§3. The viscous modes. 


Following the discussion of the solution properties of the principally inviscid modes 
considered in the previous section we now examine the effect of crossflow on the stability 
of the viscous modes whose scalings are suggested by the asymptotic results (2.15). The 
wavenumber of these vortices is predicted to be O (g»), i-e. is that found by DHS for 
the most unstable vortex in two-dimensional flow. For no crossflow, DHS showed that 
the vortex is then confined to an 0 (G~i) thick region which lies immediately adjacent 
to the cylinder y = 0 and that the spatial growth rate of the instability is 0(G* ) in the 
streamwise direction. We have, from (2.155), that when the scaled crossflow A* = O (g*) 
(where A* is defined in (2.2a)) it is to be expected that a significant deviation in the 
stability characteristics of the flow from those which occur in the zero crossflow case. 
Then, following DHS, we write 

a = k$G * , (3.1a) 

where fc 0 is 0(1) and seek disturbances confined to the layer in which 1> = 0(1) where 

$ — koG‘Y. (3.15) 

Then we consider modes for which the perturbation velocities (U, V, W) and pressure P 
take the forms 
(U,V,W,P) = 

([ 7 0 + a- -!7i+ ...,Giv 0 + GiVi + ...,G‘W 0 + Wi + ...,d Po + G-Pi + •••) 

x exp \atJ X («*) + G- *MX) + . . .) iX - iG? / (flo(<) + G- + • • •) * 

(3.2a) 

and we fix the crossflow parameter at 0(G» ) so that 

A* = AG*. ( 3 - 2& ) 

Here the unknowns U 0 , U\, V 0 , . . . , Po, Pi, • • • are functions of $ alone. 

We substitute (3.1) and (3.2) into (2.3) and compare like powers of the large parameter 
G. In this process, we need to expand the basic flow about T = 0 and so, as in (2.13), we 
assume that 

U = mi (X)Y + \ni2 (X)Y 2 + ±1X13 (X)Y 3 + . . . , (3.3a) 
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w = H2l(X)Y + \»22{X)Y 2 + lfx 23 (X)Y 3 + . 

z o 

C = 0(y 2 ), for Y < 1. 
i.From leading orders in (2.3a) we have 

/?0^0 + = 0, 

PlUo + PiUq + + ikoWi = 0 , 


and similar manipulation for (2.36) yields 

£om 11 

iX 


_ , , i«o _ Punv/^ 

<^ 2 ~ 3 k$ 2fcg/in 


dip 

+ *Ap21 = 0, 

(A<llh22 - H2lVl2)i> 2 


U ° = ^ 
Kq 


(3.36) 

(3.3c) 

(3.4a) 

(3.46) 

(3.5a) 

(3.56) 


and 


d 2 _ j /?i/ni^ 


i A 


dip 2 


fc 0 


iX 

K 0 




rr t/ 

tfi - -p-Vi 

/c 0 


| *A % ,3 , 0un2ip 2 @2Hll ip iftl 

[^(«.,23-M 21 , 13 )0 + -^ r+ ^ r _- rj 


TT jl .. 

u ° + -JT Vo ■ 


(3.5c) 


On considering the spanwise momentum equation (2.3d), at leading orders we just 
retrieve a trivial combination of (3.4a), (3.5a) and (3.56). The third order terms imply 
that 

.O . _ . 1 

A^21 


d 2 _ l + *^0 _ PlUli i> 


L dip 


c 0 




W-l - 


T 7 

Kq 


-Vi 


i A 


L 6^1lfco 


(«l«3-/.21«3)V- 3 + /,1M1 ^ 2 ' ,!!l 


2^o + 


/Cq 




m + iF ° 

Kq 


ko 

( 3 . 6 ) 

Eliminating Pi and IT i between (3.46), (3.5c) and (3.6) yields an equation for the leading 
pressure term Po of the form 

p _ , d 3 Vo dVo ( iXip Pvm\ 

Po - ko l w ~ °i? + V Ufa ( '‘ u, ‘ 22 - " 21, ‘ 12) + v ° 

, tOodVo PimvP dVo iX , % 0 dVn 

(3.7a) 
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Finally, leading order terms in the remaining momentum equation (2.3c) give that 


dP 0 

dtp 


= ko 


.dip' 


1 + 


iflo 

T 2 " 

Kq 


Pvm 


— (AU1M22 - M21M12) '* P 2 

2fcg//n 


Vo- 


XMllV> TT 

Kq 


(3-7 6) 


so that on eliminating Po between (3.7a, 6) and simplifying we obtain 


d 2 _ iClo _ Pivnip 
diP 2 ~ + If 


i A 

2fcoMll 


r 

U 11 M 22 - /U 21 M 12 ) 


(£-0 


Vo 


+ 


fcoMll 


(P11M22 - ^21M12) Vo = 


XMllV> rr 

~ P~^°- 
/c 0 


(3.8) 


It is of interest here to note that if we had chosen the crossflow to have been of the 
more obviously significant size 0(G» ) then the quadratic terms in ip in the above equations 
would be replaced by linear ones. Then /?i can be redefined to remove these terms so that 
such a weak crossflow does not have a significant efect on the viscous modes. To summarise 
thus far, the leading order vortex quantities Uo, Vo satisfy the sixth order system (3.56), 
(3.8). To simplify this pair we introduce the parameters k, /?, U and V defined by 


ko = (xdii) 7 k, /3i = X‘VhP, U 0 = nilU and Vq = X 1 V-{i v ■ ( 3 - 9a d ) 


Further, if we write 


jl — — 

flo = X'Mll^ and 


A = 


A (fillfJ-22 ~ M12M2l) 
— » 

2^11 X‘ 


(3.9 ej) 


and substitute (3.9) in (3.56) and (3.8), we obtain our final equations to determine U and 
V, i.e. 


( 


dip 2 


id (3ip i\ip 2 \ f 

1 + P “ W~ ~W ) V 

f d 2 id 

V# 2 ~ 1+ p - 


J— 

dip 2 
} c 3 


- 1 ) 9 + W V = ~W’ (3 - 10a) 

- 0 = p- (3.10b) 


\Ve need to impose suitable boundary conditions on the solutions of (3.10). First, in 
order to ensure that the disturbance is confined to the thin 0(G~ • ) zone close to the wall 
we demand that the perturbation quantities decay exponentially as ip — > 00 . Additionally, 
to satisfy the requirement of zero disturbance velocities on the wall ip = 0, U — V = 0 
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on ip = 0 and the continuity equation (3.46) further implies that dV j dip = 0 on ip = 0. 


Hence, in all, our six boundary conditions 

are 


r-r dV 



U = V = : 

dip 

cT 

II 

-S- 

a 

o' 

ii 

(3.11a) 

___ dV 



1 

l-s- 

1 

II 

u 

II 

b 

-> 0 as %l) — > oo. 

(3.116) 


Now (3.10) and (3.11) constitute an eigenproblem to be solved for the scaled vortex fre- 
quency Cl, scaled spatial growth rate j3 and scaled crossflow A. We note that once the 
eigenfunctions U , V are determined we can retrieve the spanwise velocity component of 
the vortex by applying (3.4a) and (3.9). 

The equations solved by DHS are precisely (3.10) with Cl = A = 0. They demonstrated 
that for all k we then have j3 > 0 so that the corresponding vortex flow is unstable. Further, 
as k — > 0, f3 ~ k* and as k — ► oo, (3 ~ k~ 2 . The behaviour as k — ► 0 matches with 
(2.5c) which determines the growth rate of 0(1) wavenumber vortices in two-dimensional 
flow. Since /? — * 0 as both k — ► 0 and as k — ► oo there is a most unstable vortex at 
some intermediate value of k which DHS calculated to be k = 0.476 with corresponding 
P = 0.312. The dependence of ft upon k as found by DHS is illustrated in Figure (3.1). 
Our interest is with non-zero Cl and A, but before we consider any specific calculations 
we make one observation concerning (3.10). We see that if (Cl, ft, A) is an eigenvalue set of 

(3.10) &: (3.11) then so is (—Cl* ,ft* , —A) for real crossflows A and where an asterisk denotes 
complex conjugate. Evidently we can therefore restrict our attention to postive crossflow 
parameters and do so for the remainder of the paper. As a first detailed examination of 

(3.10) , (3.11) we will study the effect of introducing the crossflow A whilst maintaining a 
zero frequency Cl — this case is treated in the forthcoming section. 


§4. Solution of (3.10) for the zero frequency case. 

We first consider the solution of governing equations (3.10) with associated boundary 
conditions (3.11) for the steady vortex problem (Cl = 0). Naturally, this problem had to 
be tackled numerically and the solution technique used in practice was based on a method 
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suggested by Malik, Chuang & Hussaini (1982). In the zero frequency work described in 
this section we solved 


M ( d2 ^ y , ^Aj? - . 

P fc3 )\j$S V + fc3 

ipU 

"■p"’ 

(4.1a) 

' d 2 P^P i\ip 2 \ p V 

<diP 2 k 3 k 3 ' P’ 


(4.1J-) 


by prescribing values for the real scaled wavenumber k and the crossflow A. Then (4.1) 
with its associated boundary conditions (3.11) were solved for the complex-valued spatial 
growth rate p. Recalling the definitions (3.2a), (3.96) we can see that when Re(0) > 0 the 
vortex is unstable and, conversely, for Re(P) < 0 the disturbance is stable. As remarked 
already we can, without loss of generality, restrict ourselves to considering A > 0. 

Malik et. al. (1982) described a boundary value method for solving eigenvalue prob- 
lems. In these methods the differential equations to be solved are reduced to a set of linear 
algebraic equations using either a finite difference discretisation or a spectral representation 
and the eigenvalues are found by solving the characteristic determinant of a generalised 
eigenvalue problem. The particular method devised by Malik et. al. uses a fourth order 
accurate (Euler-Maclaurin) finite difference scheme with nodal points distributed so as 
to resolve any singular layers. These authors used their scheme to examine the temporal 
and spatial stability of a three-dimensional compressible boundary layer flow over a swept 
wing. For further detials of this numerical method the reader is referred to Malik et. al. 

( 1982 ). 

To implement the solution procedure for (4.1) we decided to perform a sequence of 
r un s with a gradually increasing crossflow A and a fixed wavenumber k = This choice 
of k is not significant, and we duly performed further tests at other values of k in order 
to ensure that the results were qualitively similar to those we describe here. For zero 
crossflow we were, of course, just repeating the calculation of DHS presented in Figure 
(3.1) for which (3 is real and the vortex is unstable for all k > 0. For non-zero crossflows 
the spatial dependence parameter P becomes truly complex and some solutions of (4.1) are 
shown in Figure (4.1). Here we have fixed k = \ and we observe that as A increases Re(p) 
decreases so that the vortex becomes progressively stabilised. When A = 0, 0 = 0.313 
(from DHS) and Re{p) vanishes when A » 0.406, at which point Im(p) « -0.91 (see 
Figure (4.16)). As the crossflow increases yet further Re((3) becomes negative and hence 
we have the situation of a stable vortex flow. The form of behaviour of j3 when k = \ is 
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typical of that for other wavenumbers and so we may conclude, at least for these stationary 
vortices, that the effect of increasing the crossflow is to stabilise the flow. 

We then extended our computations by fixing the crossflow parameter A and solving 
(4.1) for a selection of wavenumbers k. The results obtained are illustrated in representative 
Figures (4.2a, 6). Here we plot the dependence of Re{0) upon k. In Figure (4.2a) we have 
fixed A = 0.414, i.e. a little larger than the value at which the vortex of wavenumber 
* = °’ 5 becomes neutrally stable. It is clearly seen that although Re(0) < 0 for k = 0.5, 
there is a range of wavenumbers, approximately 0.37 < k < 0.44, over which Re(fi) is stiU 
positive. Consequently, at this crossflow value there remains an interval of wavelengths for 
which the stationary vortex is unstable. However, from Figure (4.26) it is plain that when 
A is increased further (here A = 0.434) we have Re(0) < 0 for all wavenumbers and hence 
when the crossflow becomes this large the flow is completely stabilised. More detailed 

calculations showed that this occurs whenever the crossflow is greater than approximately 
0.416. 

There are two natural asymptotic regimes to investigate within the context of the 
present problem. The first deals with the large wavenumber limit k — ► oo. We do not 
address this limit here because it also arises when making large wavenumber studies of the 
temporal problem dealt with in section 5. Consequently, we shall focus our attention on 
^ * 00 tbere - Second is the case of small scaled wavenumber Jc — ► 0. We can develop 

a formal asymptotic solution of (4.1), (3.11) in this limit but rather than concentrate on 
small wavenumbers for stationary vortices alone we shall consider a slightly wider class 
of flows and shall investigate small wavelength vortices with small frequencies. Of course 
we can recover the asymptotic solution for stationary vortices without formal difficulty 
but also, much more importantly, considering small wavenumber vortices of the type just 

described will allow us to tie this analysis in with results for the temporal case in section 

5. 


§^.l Solutions of (3.10) for k «C 1. 

Guided by the eigensolutions of the problem (4.1) for small k it is found that for 

k < 1, n = 0(B) the solution vortex of (3.10) is neutrally stable for a scaled crossflow of 
size O (k-r). To demonstrate this we suppose that 

= + A = C 0 + ..., 0 = iD 0 + 
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(4.2) 



Here we have assumed that /? is purely imaginary so that D is real and with the scalings as 
in (4.2) the solution structure divides into a number of distinct regions. For the moment 
we take Cl > 0 as we will subsequently find that a modified asymptotic structure is required 

for the case Cl < 0. 

Close to = 0 we write 

rp = 09, V = V 0 +..., V = l-*Uo+.... (4-3) 

Substitution of (4.2) & (4.3) into governing equations (3.10) shows that 

( J?_ _ iCi - iD e - ice 2 ) ^ + 2iCVo = -Wo, (4.4a) 

( _ iCi _ we - ice 2 ) Uo = v 0 . (4.46) 

\de 2 ' 

We expect the disturbance to be largely confined within this O (k*) layer and so 
we demand that as e — * oo the solution decays. Inspection of (4.4) then imposes the 

conditions 

V 0 oc i + . . . , Uo OCJS+'..- M #— * “O' < 4 ' 5o > 

Additionally, boundary conditions (3.11a) demand that 

Un = V 0 = — = 0 at 0 = 0. (4.56) 

0 0 de 

Now (4.4), (4.5) together form an eigenvalue problem for the 0(1) real constants C and D. 
Before considering this further, we need to check that the disturbance decays correctly so 
as to ensure that the boundary conditions (3.116) applicable as — > oo are indeed met. 

Inspection of (3.10) subject to the scalings (4.2) suggests that when rf> = 0(1) the 
velocity components of the vortex behave according to 

V = k^U + ..., V = n*V + ..,, (4- 6 ) 

using (4.3) and (4.5a). Substitution of (4.6) into (3.10a) yields the equations 

- l) V + 2V = 0, (4.7a) 

\ dtp* ' 


and 


-iCrp 2 U = V. 


(4.76) 
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To match with (4.5a) requires V a £ as ip — > 0 and the solution of (4.7a) with this 
property is 

V = Eip* Ki(ip), ( 4 . 8 ) 

where E is some constant and K v {ip) is the modified Bessel function of order v. Solution 
(4.8) has the required exponential decay as 0 — ► oo and so, therefore, U also decays 
exponentially (see (4.76)). Hence we have obtained the complete asymptotic description 
of the solutions of (3.10) k (3.11) for neutral modes when k <C 1, Cl = 0(k$), Cl > 0. 

Technically, to complete our description of this limit the solution of (4.4), (4.5) is 
desired for various values of Cl. We have not done this here for this aspect is not within 
our principal objectives. The asymptotic work above demonstrates that for k < 1 only 
a small scaled crossflow, O (k±) is required to completely stabilise a vortex of positive 
frequency O (k*). This is an exceptionally tiny frequency and so in practice the scaled 
frequency parameter Cl defined in (4.2) is likely to be large. It can be shown (and details 
are available from the authors) that for Cl > 1 we have C = 0(51 7) and D = 0(Cl f ) and 

this information matches with further asymptotic limits to be discussed in the following 
section. 

The above describes the asymptotic solution of (3.10) & (3.11) for k < 1 and Cl smaU 
but positive. We can show however that for Cl negative a completely different structure 
must come into play, and in order to illustrate this revised case we will work with stationary 
vortices (fi=0) although our analysis may be easily modified to account for |fi| <C 1, Cl < 0. 
In the stationary case we need to solve (4.1) for k < 1. We write 

P = i0$ + . . . , A = fc»A + ..., 

(4.9) 

V = V 0 + ..., U = k>U 0 +..., 

in the region where = 0(1). Substituting these expansions in (4.1) yields 

-* + ^ 2 ) (^2 - 0 F ° + 2iXVo = (4.10a) 

+ H 2 ) Uo = Vo, (4.106) 

so that 

(w + Fo ) - 2A + ^ 2 ) Vo = 0T O , (4.10c) 
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which is an eigenproblem for the quantities p and A when we impose the inviscid condition 

V 0 > 0 as ^ — > 0 and demand that Vo — + 0 as — ♦ °o. For neutral modes we 

need P real and from the numerical work illustrated in Figures (4.1) & (4.2) we saw that 
for positive crossflows A in (4.1) the neutrally stable vortices satisfying (4.1) typically had 
Im(j3) < 0. Consequently, in terms of (4.9) we anticipate that the eigenvalues p, A of 
(4.10c) are such that p < 0, A > 0 and so (4.10c) has a critical layer structure surrounding 
^ = -pi A. Of course, the critical layer problem (4.10c) is closely related to (2.146) (in fact 
it is just a scaled version of (2.146)) which governs the large wavenumber inviscid modes 
studied in section 2. As discussed there, in order that a numerical solution on (4.10c) be 
possible a detailed analysis of the critical layer structure is required and which we will 
consider in a future article. However, the pertinent result from this critical layer problem 
is that it represents the description of how the neutrally stable vortex modes are affected 
as we pass from the parameter regime of section 2 to those considered in this section. 

Although the above description is specific to stationary vortices it is easily generalised 
to account for disturbances of frequencies ft = 00 ) and negative. The effect of the non- 
zero frequency is to introduce two extra constants into the critical layer equation (4.10c) 
but the overall properties of the problem are virtually unaffected. Again inclusion of these 
non-zero frequencies provides scope for further discussion of the critical layer problem. 

To summarise thus far, we have shown that for small frequencies the effect of increasing 
crossflow in the fundamental equations (3.10) with boundary conditions (3.11) is to tend 
to stabilise the vortex flow. Interestingly, we have also demonstrated that for small non- 
dimensional wavenumbers the structure of the neutrally stable modes is dependent largely 
on the sign of the frequency of the modes. For frequencies ft > 0 the vortex structure is 
dictated by the solution of eigenproblem (4.4) whereas for ft < 0 we obtain a critical layer 
structure which matches with the large wavenumber inviscid modes described by (2.14). 
Having studied the low frequency problem we now move on to consider larger frequencies, 
indeed we examine the case for which the scaled frequency ft = 0(1). 
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§5. The vortex problem (3.10) for 0(1) frequencies. 

Here we examine the properties of Gortler vortices of wavenumber 0(G») and fre- 
quency 0(G* ) in the presence of crossflow of size 0(G» ). In this case, of course, the 
scaled parameters Cl, 0 and A in (3.10) are all order one quantities. For the purpose of this 
section we will concentrate almost exclusively on neutrally stable states for the results of 
the previous sections have suggested how the stability of the flow will be altered by chang- 
ing particular parameter values which give rise to neutral modes. To solve for neutral 
modes we again used the method described in section 4 subject to a few slight changes. 
We then specified the 0(1) frequency Cl and the wavenumber k in (3.10) and treated these 
coupled ordinary differential equations as an eigenproblem for the crossflow parameter A 
and the spatial development quantity (3 which was supposed to be purely imaginary. In 
this way we ensured that the resulting eigenvalues corresponded to neutrally stable modes. 

Some sample results are illustrated in Figures (5.1) - (5.3) where we show the depen- 
dence of the crossflow A and the parameter 0 on the scaled wavenumber k for the three 
frequencies Cl — 1 , 0 and — 1 . Briefly, we can see that neutral modes appear to be possible 
over the complete range of wavenumbers k and that for large k the crossflow needed to 
produce neutrally stable modes is quite small. However, a striking difference between the 
neutral modes is observed depending on whether Cl > 0 or Cl < 0. In the former case (see 
Figure (5.1a)), as k — ► 0 A tends to an 0(1) value whereas in the latter case A — ► 0 as 
k — ► 0 (Figures (5.2a, 5.3a)). A discussion of the asymptotic properties of these two cases 
is given later. We note here that on recalling the results of section 4 together with the 
zero crossflow calculation of DHS we can infer that for crossflows greater than the neutral 
mode values on Figures (5.1a — 5.3a) the corresponding vortex motions are stable whilst 
the opposite is true for lesser crossflows. 

We can be confident that the results shown in Figures (5. 1-5.3) are sufficient to demon- 
strate that increasing crossflow again stabilises these non-stationary vortex modes in a 
similar way as was described for the stationary case in section 4. As before, there exist 
several regimes in which asymptotic description of (3.10) & (3.11) is possible. We first 
concentrate on properties of this eigenproblem in the high wavenumber limit, k » 1. 
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§5.1 The high-wavenumber ( k 1) limit 

In this limit we are concerned with obtaining an asymptotic solution of the system 

( 3 . 10 ) , ( 3 . 11 ) for fc > 1, H = 0(1). We can initially consider the problem in the context 
of neutrally stable modes in the spirit of the calculations just described. It has been seen 
that for low wavenumber modes the vortex structure appears to be critically dependent 
on the sign of the frequency Cl but inspection of Figures (5.1-5.3) suggest that this is 
unlikely to be the case for high wavenumber modes and consequently there is no need to 
specify Cl further than taking it to be of order unity. An idea of the expected structure 
can be gleaned by studying Figure (5.4). Here we present the eigenfunctions of (3.10) & 

(3.11) when this pair of equations is solved for neutral modes with Cl = 0 and k « 1.35. 
Unfortunately, for values of k much larger than this value we found that it was difficult 
to obtain satisfactory numerical convergence but the evidence of Figure (5.4) and close 
inspection of the governing equations are sufficient in order to strongly suggest that the 
disturbance is confined to an 0(1) thick region which moves away from 4> = 0 as k — > oo. 
If we write 

Tp = ck 5 +'4>, (5.1a) 

p = i {poP + (hk- 2 + . . .) , ( 5 - 16 ) 

and 

A = Xok~ 2 + Aifc“ 7 + ..., ( 5 - lc ) 

together with the eigenfunction expansions 

U = Uo + k~ 5 Ui + . .. , v = k 2 Vo + fc _3 Vi + ..., ( 5 - lc 0 

and substitute (5.1) into (3.10) we find the following. At leading orders we obtain 

00 = — Ao c, ( 5 - 2 ) 

and, at next order, 

- 1 + i0 - i (Pi c + Aic 2 )] \—j 2 - l] v 0 = -cUo, 

L d0 2 J L ^ 2 (5.3) 

— 1 + i0o4> — i (P\c + Aic 2 ) I Uo = Vo. 

L # 2 
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By redefining tp — ip — (( p\c + Aic 2 )//?o) we simplify these equations to 


- 1 + 10001 I -ix - 1 

dip^ dip^ 


V 0 = -cU 0 , 


L dtp 5 


- 1 + iPoi> 


U 0 = Vb, 


which need to be solved for real Pq and c subject to the boundary conditions 

dV o 


V o, Vo, 


dip 


0, as ip — ► ±oo. 


(5.4a) 

(5.46) 

(5.5) 


The precise solutions of (5.4) &: (5.5) are relatively unimportant- the crucial conclusion 
we can draw from this asymptotic work is that for an 0(1) frequency we only need a very 
small 0(fc” 2 ) crossflow to stabilise modes with wavenumber k 1 . In passing, we notice 
that the eigenproblem (5.4) may be used to deduce other results. In particular, we have 
the result that the description of the asymptotic problem k — > oo is independent of the 
order one value of f 2 a result which is in keeping with the numerical evidence presented 
earlier. Indeed, the above is valid for all frequencies Fl = o(k 2 ) so that in particular it is 
valid for the high wavenumber limit of the zero frequency work examined in section 4. 

We now make some comments about the low wavenumber limit, k — ► 0. 


§5.2 The low wavenumber (k 1) limit for Fl > 0 


We have already observed that different eigenproperties of (3.10), (3.11) are expected 
for fc < 1 and that the sign of Fl is crucial in determining these properties. Consequently, 
here we take Fl > 0 and examine the case Ft < 0 later (the zero frequency case has already 
been dealt with). When Fl > 0, k <C 1 the eigenfunctions of (3.10) are concentrated in a 
thin region near ip — 0, see Figure (5.5). In this figure we illustrate the eigenfunctions V , 
V for the case Fl = 1, k = 1.3 x 10~ 2 . This suggests seeking a multi-zoned asymptotic 
solution structure of (3.10) for k <C 1, Fl = 0(1), and we now show how this may be 
achieved. Initially, we suppose that the crossflow A and spatial dependence j3 assume the 
forms 

A = Ao + X\k 4 -f A 2 k 3 + . . . , 

(3 = i \(3ok* + pi k* + p2k + . . .] , 
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(5.6a) 

(5.66) 



and we recall that we are seeking neutral disturbances so that the 0(1) constants A 0 , Ai, . . . , 
00, ft,... are real. After experimentation, we can infer that the disturbance is largely 
trapped in an 0(S*) thick region which is located at a distance O (**) from iP = 0. 


Guided by this, when 

ip = k* ipo + k* ip, 

(5.7a) 

say, with ipo — 0(1), we need to expand 


U = U 0 k~ » + Uik~ « + U 2 + • • • , V = Vo + Vifcr + V 2 k » + . 

(5.76) 

Substituting (5.6), (5.7) in (3.10) yields that in this region (refer to it as 

zone I), 

ft + Po^Po + AoV>o = 

(5.8a) 

00 + 2ipoXo = Mo + Mo = °- 

(5.86, c) 

If we define ft = Mo + A 2 V>o " (£l/ 4A o) and write 




we deduce that within zone I the leading order disturbance quantities Uo, V 0 satisfy 


— -j + 2iAoVo = -ipoUo, (5.9a) 

dip 

[ - iCl - i\olp 2 } ft) = Vo. ( 4 5 - 96 ) 

l dip J 

For the disturbance to be effectively confined to zone I we require the solutions to 
decay as ip * ±oo. We can show that the potentially most dangerous mode is even in ip, 

and in addition it is easy to prove that 

ft) * — T3> as 1^1 ><x> * (5 ' 10) 

ip Ao rp 

la the course of his investigation into Gortler vortices at 0 (G») wavenumbers in 
three-dimensional boundary layers, Hall (1985) considered the eigenproblem (5.9), (5.10). 
He found that the relevant eigenvalues are 

4 = 4.71 and A = -3.59. (5.11a, 6) 

Aq a o 


— iCl — iXoip 


A.U 
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Even though we now have enough information to determine Ao and via (5.8a, 6) and 
(5.11a), we must check that this zone I solution obtained above can be satisfactorily 
matched to the requisite boundary conditions (3.11). Since zone I is fixed at a distance 
0 {k 3 ) from the wall ip = 0 it is natural to define a second region, zone II , by 


ip = k*9. 


(5.12) 


By (5.7 6) and (5.10), within II the eigenfunctions must expand according to 


V - vok* + . . . , U = uok* + (5.13) 

Substitution of (5.6), (5.12) & (5.13) into the governing equations (3.10) yields, on appli- 
cation of the results (5.8), that 


(9 - ip 0 ) 



— 2vq — 0, 


(5.14a) 


and 


2 

iAo (Q ~ ipo ) uq = -u 0 . 


(5.146) 


Strictly, of course, we need to consider the regions 0 < ip 0 and 9 > ip 0 separately, and we 
shall concentrate on the former first. 


The region 6 < 0g 


We first recall that the boundary conditions (3.11a) need to be satisfied at 9 = 0 
and, guided by result (5.10), we solve (5.14a) subject to the conditions vq — > j + . . . as 
9 — ► vo = 0 at 6 = 0. This yields 


fp - (jpQ - fl ) 3 

V>o Wo ~ 0) 


(5.15a) 


and, by (5.146), 

* ko - (^o - 0) 3 ] 

(5 - 156) 

We see that (5.15) satisfies only two of u 0 = v 0 = = 0 at 9 = 0 and to force the third 

condition to be met we now add a viscous layer, zone III, at 9 = 0. If in this layer we 
define the 0(1) co-ordinate 


<P = k(, 


(5.16) 
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and the expansions 


U = uok * + ■ • • > 


V = hok* + 


(5.17) 


then to satisfy (3.11a) and to match with (5.15) we require 


* a U uil A + n 

U 0 - vo = -JT = 0 at £ — 0 

at, 


(5.18a) 


A. V 

4o “ * Wif + 


> "75 + * • * ? as £ > 00 ’ 

00 


Substituting (5.16), (5.17) into (3.10) yields 


(|*- in )§? = 0 ' (^- m ) 4 » = 40 ' 


The solutions of (5.19) satisfying (5.18) are 




(5.186) 


(5.19) 


4 ° = [* - exp ( _(1 + iKv/ ?) + Sv 4 + m 

where m 2 = tft and JJ is chosen so that fi 0 (0) = 0. In this way we have shown that 
the solution proposed so far satisfies the required conditions at 0 = 0. We now need to 
reconsider the 0(S*) zone, zone II, where 0 = kU and 9 > 0o- Here expansions (5.13) 
and equations (5.14) again apply but this time we need to solve (5.14) so that n 0 ♦ j + • • • 
as $ — ► and uo — * 0 as 9 — * oo. The required solution is simply 


(9 -0o)’ 


(5.20a) 


and. then 


„ 0 1 3 - (5 ' 20,>) 

Ao (9 - 0o) 

To complete the structure we are led to consider a fourth zone, IV, where 0 = 0(1) 
and from (5.20) this means that the eigenfunctions have behaviours of the forms 


V = *o k* + . • • , U = uok * + 
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To match with (5.20) requires wo — * j + . . . and u 0 — * + ... as ip — ► 0. Substitution 

into (3.10) finally gives that 

O d WO n 

^ - ( 2 + ^ 2 )®0 = 0, (5.22a) 


UQ = 


»'f>0 

Aq0 2 ' 


(5.226) 


The solution of (5.22a) with the required algebraic growth as ip — > 0 together with 
exponential decay as ip — > oo is 


voocip>Ki(ip), (5.23) 

where K v {ip) is the modified Bessel function of order u. Clearly, &o also has exponential 
decay as ip — > oo and so we have a complete asymptotic solution of (3.10), (3.11) for 
^ — 0(1), k < Cl and p imaginary. The full four zoned structure of this solution is 
summarised in Figure (5.6). Returning to (5.8a, b) and (5.11) yields that 

Ao~0.46f2«, ~ 1.4717* and /?o « — 1.3617* . (5-24) 

Hence we conclude that for small wavenumber vortices, a crossflow of approximately 
0.4617* + ... is required to stabilise the flow. We note also that as — ► 0 in (5.24) 

our predictions match with the solution of the k 1, Cl — 0 (^ 7 ) problem described in 
section 4.1. Having now obtained the solution structure for k <C 1, Cl > 0 we turn now to 
consider the second case, Cl < 0. 


§5.5 The low wavenumber (k <C 1 ) limit for 17 < 0 

Having deduced the difference in behaviours of solution properties of (3.10), (3.11) 
when k < 1 depending on the sign of Cl we now study the case 17 < 0. From Figure 
(5.2) we conclude that the most fundamental difference between the two cases is that for 
17 < 0 both p and A tend to zero as k — > 0 in order that neutrally stable modes may be 
obtained. In order to study the structure here we show the eigenfunctions of the problem 
(3.10) when 17 = -1 and k = 8.7 x 10~ 3 , see Figure (5.7). Clearly the eigenfunctions 
are confined to a thin region centred away from the wall ip = 0 and suggests that some 
critical layer type structure is operational. This is in agreement with the low frequency 
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work of the previous section where we proved that when 0 is negative a similar critical 

layer approach is required, see equations (4.10). 

To describe the detailed structure we suppose that across the majority of the flow, 

where xp — 0(1), we write 

V = V 0 +..., U = k>U 0 + .... (5.25a) 

and suppose that the streamwise dependence parameter 0 and scaled crossflow A take the 
forms 

0 = i (fiok* + . . .) and A = (Ao k* + • - •) • (5.256) 

Substitution in (3.10) yields that Vo satisfies 

(poxp + AoV> 2 ) ~ Vo) ~ 2 ^0 {&o4> + A^ 2 ) Vo = (5.26) 

This equation, which is identical to (4.10c), needs to be solved subject to V 0 = 0 at xp = 0 

and V Q > o as xp — ♦ oo. As before, there is a critical layer structure at xp = -0o/ A 0 which 

needs further study before a full numerical solution of (5.26) can be effected. However, we 
can draw the general conclusion that for H < 0 the crossflow needed to maintain neutrally 
stable vortices vanishes as k — > 0. 

In addition to the work on neutral modes described here we have performed a few 
calculations for non-neutral modes. In Figure (5.8) we show the effect of crossflow on 
the stability of the vortices. We considered the two frequencies 0 = -1,1 and took the 
crossflow parameter A = ,1,.2,.3,.4 and, for various wavenumbers k, solved (3.10), (3.11) 
for the complex valued 0. Figure (5.8) shows the growth rate curves for vortices with the 
above crossflow parameters . We can see (as has already been commented upon) that for a 
fixed wavelength vortex an increasing crossflow progressively stabilises the flow. Further- 
more we see that instability occurs only for a finite range of values of the wavenumber; 
and we note that the width of this unstable band decreases as the crossflow increases. We 
also note that the flow is more unstable when the imposed frequency is positive. We now 
conclude with some further observations and some discussion. 
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§6. Conclusions. 

In the previous sections we have described the process by which the Gortler vortex 
instability mechanism is destroyed by the introduction of a weak crossflow into a centrifu- 
gally unstable basic state. We choose to concentrate on the situation when the crossflow 
is of size Re~ » because that is the crucial size at which the crossflow first has a significant 
effect on Gortler vortices. In addition we elect to concentrate on the large Gortler number 
limit because that limit describes inviscid Gortler vortices and nonparallel effects can be 
accounted for using asymptotic means. We found that at order one wavenumbers the in- 
troduction of a crossflow causes the growth rate-wavenumber curve of inviscid theory to 
split into two distinct parts. Thus above a critical size of the crossflow a finite band of 
unstable modes exists to the right of a = 0. This band of modes takes on the asymptotic 
structure of crossflow vortices when the crossflow is increased further. The growth rates 
of these modes increases monotonically with the crossflow parameter and will become for- 
mally of order Re * when the crossflow is comparable to the streamwise basic velocity field. 
The second branch of unstable modes occurs at progressively higher wavenumbers as the 
crossflow increases. The growth rates of these modes at any fixed value of the crossflow 
increases like a> for large a so it would appear at first sight that, at a fixed value of the 
crossflow, the viscous modes have the highest growth rates. However this is not the case 
since viscous effects, as in DHS, serve so as to stop this monotonic growth of j3 with a. In 
the absence of any crossflow the growth rate peaks at wavenumber of size G» and has a 
maximum of order G* . At higher values of a the growth rate tends to zero and the right 
hand branch neutral structure of Hall (1982a) is retrieved. In the presence of a crossflow 
we showed that the second, semi-infinite band of unstable modes predicted on inviscid 
grounds becomes reduced to a finite band at high wavenumbers when the crossflow is for- 
mally of order G« . At this stage the growth rates of the crossflow vortices associated with 
the modes to the right of a = 0 have growth rates of size G» comparable with those of 
the viscous modes. However as the crossflow is increased further the growth rates of the 
crossflow vortices increase whilst the viscous modes become stabilized. We further note 
that the unstable band of modes which emerge to the right of 02 connect with the viscous 
modes in an inviscid manner; by this we simply mean that their structure is governed by 
a critical layer scenario near a?. 

In Figure (6.1) we have sketched the growth rate-wavenumber dependence determined 
above for the three situations: (a) A* (b) A* ~ G« , (c) A* > G» . In (a) the crossflow 
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is strong enough only to alter the in viscid modes which have growth rate of size ; the 
viscous modes have growth rates of order G» and are not affected by the crossflow. In (b) 
both modes have growth rates of order but the inviscid modes are now independent of 
the wall curvature at zeroth order; thus the inviscid modes have deformed into crossflow 
vortices. The viscous modes at this stage depend on the crossflow and become completely 
stabilized at a finite value of A * /G*. However before this occurs there are two points 
at which the growth rate has a local maximum. The relative importance of the maxima 
associated with the viscous and inviscid modes switches over as the crossflow increases. In 
(c) only the crossflow vortex mode remains and at zeroth order it is independent of the 
curvature of the wall. 

The question of which of the possible modes discussed above is the most likely to be 
observed in an experiment can only be answered by a receptivity calculation of the type 
given by DHS for Gortler vortices in two-dimensional boundary layers. In fact the analysis 
given in Secion 6 of the latter paper is easily modified to discuss the receptivity problem 
for the viscous modes described by (3.10). Since this type of disturbance is concentrated 
near the wall then, following the analysis of DHS, it is a relatively simple matter to show 
that the ‘coupling’ coefficient associated with a roughness element at the wall is 0(1). By 
this we simply mean that the streamwise velocity component of a mode stimulated by an 
obstacle at the wall is comparable to the fluid velocity induced near the obstacle. However, 
as in DHS, the coupling coefficient depends on the vortex wavenumber and, if all modes 
are excited, sufficiently far downstream of the obstacle the fastest growing linear mode will 
dominate the flowfield. We shall concentrate on the case (b) above since by taking the 
appropriate limits in that situation (a) and (c) can be recovered. 

We suppose then that the obstacle stimulates disturbances on a spanwise lengthscale 
comparable to the boundary layer lengthscale; this means that the inviscid modes satisfying 
(2.8a) wifl be stimulated. We shall now briefly describe how the coupling coefficient for 
this type of mode can be derived. 

Suppose then that the roughness element defined by 

Y = h(X,Z), (6- 1 ) 

generates a disturbance field governed by (2.3). We assume that there are no oscilla- 
tions in the free-stream so that the disturbed flow will be steady. Following DHS it is 
straightforward to show that the system (2.3) must now be solved subject to 

U = —u Y h, V = 0, W = -w Y A*, at Y = 0 (6.2) 
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instead of U - V - W - 0 at Y = 0. More precisely we note that we must replace 
ia by a partial Z derivative in (2.3) unless the X and Z dependences of h are separable. 
For simplicity we asume that this is the case so that without loss of generality we can 
take h to be function of X alone in the following discussion. However we note that, if 
we wish to generate the flowfields in a neighbourhood of the obstacle, the particular Z 
dependence of h must be accounted for when inverting the Fourier transform (in Z) of the 
disturbance. Here we are merely interested in finding the coupling coefficients associated 
with stationary vortices of wavelength comparable with boundary layer lengthscale so it is 
sufficient to take h = h(X). 

Since the forcing takes place at the wall we must in effect determine how an inviscid 
disturbance satisfying (2.8a) adjusts in a viscous wall layer so as to satisfy the correct 
boundary condition at the wall. It is easy to see that this layer must be of thickness G~ • 

so that (6.2) suggests the following expansions for (U, V,W,P) in the wall layer where 
Y = G»q, 


( U > W,P) - ( U° + . . . , G» V° + . . . , Gj W° + . . . , G> P° + . . .) (l 4- Q(G~ « )) 


(6.3) 


The X dependence of the disturbance field is now taken care of using a Laplace transform. 
We suppose that h varies on a G’i lengthscale in the X direction. Thus we write 


h = 


/ 0 , 

ltf(GT[jr-Xo]), 


if X < X 0 , 
if X > Xq , 


and determine the flow induced beyond X = X 0 . Thus we now let /? be the Laplace 
transform variable with respect to the variable G> [X - Xo]. After a little work we can 
show that in the wall layer U° and W° are related by 


t GO 

/H7° + ialV 0 = 3hN* j Ai(N^q)dq, 
where h is the transform of H and 


(6.4a) 


N — P/J-ll + 21i 

whilst Ai is the Airy function. The function V® is then given by 

F° = -3hN f ( / Ai(N>b)db)dq. 

0 , 


(6.46) 
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We see that the above solution has V° tending to a constant V° when q — » oo so that in 
the main part of the boundary layer the appropriate expansions of U, V, W are 

(17, V, W) = (G~*Uo + Vo+ • ..,G» Wo + . . .) {l + 0(G~>)) . 

The function V 0 is then found to satisfy (2.8a); but the appropriate boundary conditions 
for this equation now become 


V 0 = V^J , at Y = 0, Vo — ♦ 0 as Y > oo 


For given values of the transform variables 0 and a this inhomogeneous differential system 
for V 0 can be solved. However if 0 = 0 an eigenvalue of (2.8a, b, c ) then the solution of the 
forced problem near 0 must be found by writing 


V 0 = 


Voo ^ 

( 0 - 0 ) 


(6.5) 


where V 00 = AV 0 . Here V 0 is the eigensolution of (2.8a, 6, c) with 0 = 0. The constant A 
may be written as ^ v O^( 0) 

A = /“ M(q)V 0 (q)dq 

Here the function M is obtained by differentiating the left hand side of equation (2.8a) 
with respect to 0 and then setting V 0 = % Having determined V 0 we then see that in 
the main part of the boundary layer the disturbed downstream velocity component in 0, a 

space near 0 = 0 is given by 

_ — G~ « AVptty (6.6) 

_ (0 - 0){/3u + i\aw)' 

It follows that when the Fourier-Laplace transform U{0, a,Y) is inverted then the pole at 0 

will lead to exponentially growing eigensoluti. proportional to AG • e p ( Thus 

the coupling coefficient for inviscid vortex modes is of order G~ • so that wall roughness is 
a more efficient stimulator of the viscous modes which, from DHS, have an 0(1) coefficient. 
However as the size of the crossflow is increased the growth rates of the inviscid modes 
increases to be larger than the viscous growth rates so this slight difference in size of 
the coupling coefficients will be irrelevant beyond a critical size for the crossflow and the 
observed instability will be a Rayleigh instability rather than a centrifugal one. 
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Finally we say a few words about the relevance of the above calculations to practical 
flows. In the obvious practical situations where Gortler vortices are thought to be a likely 
cause of transition the basic state is three-dimensional. Thus for example there is little 
doubt that the flow over a turbine blade is three-dimensional and our results suggest that 
it could therefore not support Gortler vortices. On a swept wing the basic state is three- 
dimensional and, for order one values of the angle of sweep, we can be confident that Gortler 
vortices cannot exist. However though the size of crossflow required to formally destroy 
the mechanism is of size Re~ > G* which is necessarily small because of the definition of 
the Gortler number the numerical constant multiplying this factor will be a function of 
the particular basic state. It could well be that for some swept-wing flows the Gortler 
mechanism might survive to modest angles of sweep appropriate to practical situations. 
Clearly this matter should be the topic of a careful experimental investigation of the Gortler 
mechanism in three-dimensional flows. 
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Figure (2.4a,b). Imaginary parts of the amplification parameter (3 for various 
scaled crossflow A. a) A = 0,1,2, 3 and 4, b) A = 4.5,4.75 and 5. 
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Figure (2.4c). Imaginary parts of the amplification parameter $ for crossflow sizes 
A = 5, 10, 15, 20 and 25. For each crossflow the cross indicates the classical value of such 
that (2.12) is satisfied and so neutrally stable solutions of the GSW crossflow instability 
problem arise. 
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Figure (5.1). Values of crossflow A and parameter 0 needed in o.^er to ensure that 
vortices of frequency ft = 1 are neutrally stable, a) A and b) 0 as functions of vortex 
wavenumber k. 




Figure (5.1). Values of crossflow A and parameter /? needed in order to ensure that 
vortices of frequency ft = 1 are neutrally stable, a) A and b) B as functions of vortex 
wavenumber k. 
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Figure (5.2). Values of crossflow A and parameter f3 needed in order to ensure that 
vortices of frequency 0 = 0 are neutrally stable, a) A and b) (3 as functions of vortex 




Figure (5.2). Values of crossflow A and parameter 0 needed in order to ensure that 
vortices of frequency Cl = 0 are neutrally stable, a) A and b) 0 as functions of vortex 
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Figure (5.4). Eigenfunc 
= 1.35, n = 0. a) Re(i7) 
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Figure (5.4). Eigenfuncti 
= 1.35, = 0. a) Re(£7), 
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Figure (5.4). Eigenfuncti 
= 1.35, Cl = 0. a) Re(U), 
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Figure (5.4). Eigenfunctions U, V of (3.10) & (3.11) f( 
= 1.35, Si = 0. a) R e(U), b) Im(ff), c) Re(V r ), d) Im(V) 
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Figure (5.5). Eigenfunctions U , V of (3.10) & (3.11) for neutrally stable vortices with 
= 1.3 x 10- 2 , Cl = 1. a) R e{U), b) c) Re(tO, d) 
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Figure (5.5). Eigenfunctions U, V of (3.10) & (3.11) for neutrally stable vortices with 
= 1.3 x 10- 2 , ft = 1. a) Re(ff), b) c) Re(V), d) Im(P). 
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Figure (5.6). Schematic diagram of asymptotic structure of the solution of (3.10) & 

(3.11) for the case of small vortex wavenumber i<l. The solution configuration divides 
into four distinct zones: Zone (I) is a thin region of depth 0(k* ) at a distance 0(k j ) from 
the wall, see definition (5.7a). This zone is embedded within (II), a region of thickness 
O(fc’), see (5.12). The structure has an additional thin viscous layer at the wall, Zone 
(III), which enables boundary conditions (3.11a) to be met. Finally, an outer region (IV) 
of depth 0(1) facilitates the exponential decay of the disturbance solutions far from the 
wall, see (5.21). 
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0.10 0.15 0.20 0.25 0.30 0.35 ( 

Figure (5.7). Eigenfunctions U, V of (3.10) & (3.11) 
= 8.7 x 10- 3 , Cl m -1. a) R e(U), b) Im(U), c) Re(V) 







Figure (5.7). Eigenfunctions U, V of (3.10) & (3.11) for neutrally stable vortices with 
= 8.7 x IQ" 3 , n = -1. a) R e(U), b) Im (U), c) Re(V), d) Im(V). 
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